Uncorrelated Estimates of Dark Energy Evolution 
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Type la supernova data have recently become strong enough to enable, for the first time, con- 
straints on the time variation of the dark energy density and its equation of state. Most analyses, 
however, are using simple two or three-parameter descriptions of the dark energy evolution, since 
it is well known that allowing more degrees of freedom introduces serious degeneracies. Here we 
present a method to produce uncorrelated and nearly model-independent band power estimates of 
the equation of state of dark energy and its density as a function of redshift. We apply the method 
to recently compiled supernova data. Our results are consistent with the cosmological constant 
scenario, in agreement with other analyses that use traditional parameterizations, though we find 
marginal (2-cr) evidence for w(z) < — 1 at z < 0.2. In addition to easy interpretation, uncorrelated, 
localized band powers allow intuitive and powerful testing of the constancy of either the energy 
density or equation of state. While we have used relatively coarse redshift binning suitable for the 
current set of ~ 150 supernovae, this approach should reach its full potential in the future, when 
applied to thousands of supernovae found from ground and space, combined with complementary 
information from other cosmological probes. 



I. INTRODUCTION 

Recent measurements of the distance-redshift relation 
using type la supernovae (SNe la) 0, obtained using 
the Hubble Space Telescope further strengthened the evi- 
dence that the rate of expansion of the universe is increas- 
ing in time 3]. This accelerated expansion is ascribed 
to a mysterious component called dark energy that com- 
prises about 70% of the energy density of the universe. In 
addition to supernova data, additional pieces of evidence 
come from the combined study of the large scale struc- 
ture and the cosmic microwave background anisotropy 
measurements (If . While the presence of dark energy is 
by now well established, we are at an early stage of study- 
ing and understanding this component. It is hoped that 
more accurate cosmological measurements will further 
constrain parameters describing dark energy and even- 
tually shed light on the underlying physical mechanism. 

Dark energy is most simply described by its present 
day energy density relative to the critical value, £Ide, 
and its equation of state defined as the ratio of pressure 
to density, w = pde/pde I n general, w is allowed to 
freely vary with time (or redshift), as is pee- In practice, 
it is difficult to constrain w(z) or, say, the scaled energy 
density f(z) = Pde{z)/ Pde(0), when they are described 
by more than a few parameters due to severe parameter 
degeneracies entering the observable quantity (luminos- 
ity distance, in the case of SNe la). Even though it is in 
principle possible to recover the function f(z) or w(z) 
directly from supernova measurements [||, in practice 
one has to fit the noisy data with a smooth functional 
form |jj which introduces error and bias (for a valiant 
attempt to do this with current data, see ||). Another 
general approach is to model w or / using a cubic spline 
in redshift (e.g. @), but again the paucity of data lim- 



its the spline to a few points in redshift, while having 
more points would correlate the measurements making 
the interpretation somewhat difficult. 

Constraints from the new SN la data Q suggest that 
dark energy is consistent with the cosmological constant 
scenario QjjJ , agreeing with previous work [hJO] . How- 
ever, these (and other) analyses are typically based on 
particular models — either a linear variation with red- 
shift or the evolution that asymptotes to a constant 
w at high redsh ift [RJ , or perhaps a more complicated pa- 
rameterization [14( — that are used to describe redshift 
variation of the dark energy equation of state. While 
these forms do a very good job in fitting w(z) due to a 
variety of proposed mechanisms that could be responsible 
for dark energy |15| , one should keep in mind that we are 
far from having any solid leads as to what to expect for 
the dark energy evolution. Given the constant increase 
in the quality and quantity of SN la data, it is timely 
to consider whether one can use current data to derive 
model independent conclusions on the evolution of dark 
energy. 

In this paper, we introduce a variant of the principal 
component analysis advocated in Ref. . We make use 
of the most recent type la supernova data from Ref. 
and present a view of dark energy complementary to 
other approaches. At the same time, we are seeking to 
answer one of the most important questions at present: 
is dark energy consistent with the cosmological constant 
scenario or not? Our analysis is facilitated by the fact 
that our measurements are completely uncorrelated. Fi- 
nally, we briefly comment on the applicability of this ap- 
proach to future datasets. Throughout we assume a flat 
universe. 
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II. METHODOLOGY 

We would like to impose constraints on the parameters 
Pi (i = 1 ... N) that describe the dark energy equation of 
state w(z) or its energy density f(z), each pi being suit- 
ably defined in the i th redshift bin. In addition to these, 
we have two more parameters: matter density relative 
to the critical Qm = 1 — ^de and the Hubble constant 
h = i/o/(100km/s/Mpc). We first marginalize the full 
(N + 2)-dimensional likelihood over these two (for the 
priors and assumptions, see Sec. Illlfl . and project them 
onto the pi space. The covariance of the N resulting 
parameters is 

C^( P p T >-(p>(p T ) (1) 

where p is the vector of parameters pi and p T its trans- 
pose. These parameters can now be rotated into a basis 
where they are diagonal by choosing an orthogonal ma- 
trix W so that it diagonalizes the Fisher matrix 

F = C 1 = W T AW (2) 

where A is diagonal. It is clear that the new parameters 
qi, defined as q = Wp, are uncorrelated, for they have 
the covariance matrix A -1 . The qi are referred to as the 
principal components and the rows of W are the win- 
dow functions (or weights) that define how the principal 
components are related to the pi. We refer the reader to 
Huterer & Starkman ^| for a discussion on the applica- 
tion of principal components to the dark energy equation 
of state. 

Let us now define W by absorbing the diagonal el- 
ements of A 1 / 2 into the corresponding rows of W, so 
that W T W = F . Then, as emphasized by Hamilton 
and Tegmark |l7j in the context of matter power spec- 
trum measurements, there are infinitely many choices for 
the matrix W, as for any orthogonal matrix O, OW is 
also a valid choice that makes the parameters qi uncorre- 
lated. While the principal components, qi, have several 
nice features — in particular, the best-determined qi are 
smoother and have support at lower redshift than the 
poorly determined ones — their corresponding window 
functions are oscillatory making the intuitive interpreta- 
tion of the components somewhat difficult. 

Here we advocate another choice for the weight matrix 
W: th e sq uare root of the Fisher matrix, W = F 1 / 2 = 
C -1 / 2 17] . This choice is interesting since the weights 
(rows of W) are almost everywhere positive, with very 
small negative contributions, and this has been recog- 
nized as a useful basis in which to represent measure- 
ments of the galaxy power spectrum from large-scale 
structure surveys |18j . The matrix W is computed by 
first diagonalizing the Fisher (inverse covariance) matrix, 
F = O t AO, and then defining W = T A 1 / 2 0. We 
normalize W so that its rows, the weights for pi, sum to 
unity. With this choice, Eq. Q shows that the covariance 



of the new parameters, q = Wp, is 

(to- W) («-(»») = £(F v^ (F vV <3) 

a b 

and parameters qi are manifestly uncorrelated. Further- 
more, their weights are mostly positive and are localized 
in redshift fairly well. We illustrate this in the next sec- 
tion using current supernova data. 

III. RESULTS 

We perform the analysis of the "gold" dataset from 
Riess et al. Q. First, we need to parameterize w(z) and 
f(z) in redshift, thereby defining the parameters pi from 
Sec. |H] We choose w(z) to be piecewise constant in red- 
shift and f(z) to be piecewise linear (and continuous). 
These two assumptions are consistent, since the two func- 
tions are related as w(z) = 1/3 (l+z)f'(z)/f(z)— 1. Note 
that, in the limit of the large number of parameters pi, 
the shape of the function across the redshift bin becomes 
irrelevant. 

We choose N — 4 bins with redshifts < z < 0.2, 
0.2 < z < 0.4, 0.4 < z < 0.6, 0.6 < z < 1.8, for both w{z) 
and f(z) constraints. While our choice of the number of 
parameters (or bins) is limited by the computing power 
required to perform the maximum likelihood analysis, we 
have repeated the same analysis with five parameters in 
each case and found consistent results. Future SN la data 
will lead to better constraints at all redshifts, requiring 
more parameters and perhaps the use of Markov chain 
Monte Carlo techniques, but for our purpose a simple 
analysis is sufficient. Furthermore, we have explored in 
detail the choice of the redshift binning, trying to strike 
a balance between band powers being narrow and hav- 
ing small error bars. Not surprisingly, we find that the 
constraints on w(z) or f(z) are much better at low red- 
shift, and we put three of our four bins there, choosing 
their widths so as to get comparable constraints in each. 
We have varied the exact spacing of the bins, and found 
results consistent with the same underlying w(z) or f(z). 

Finally, we describe the piecewise linear f(z) as follows: 
we write f(z) = 1 + g(z) (note that g(z) = corresponds 
to the cosmological constant scenario). We describe g(z) 
by the sawtooth basis in redshift, where each tooth is 0.2 
wide and peaks in the middle of the corresponding bin. 
The highest-redshift bin presents a problem, since it is 
much wider than the others and implies that f(z) may 
be forced to vary strongly across this bin. To prevent this, 
we make all basis vectors of the sawtooth 100% correlated 
in the highest redshift bin, essentially making f(z) flat 
across this bin. We have checked that these details do not 
affect the results appreciably by repeating the analysis 
with a few alternative choices, and we believe that these 
assumptions are reasonable and intuitive. 

The analysis is now straightforward: we compute the 
goodness-of-fit statistic \ 2 f° r each model in the six- 
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FIG. 1: Uncorrelated band-power estimates of the equation of state w(z) of dark energy are shown in panel (a). Vertical 
error bars show the 1 and 2-a error bars (the full likelihoods are shown in panel (c)), while the horizontal error bars represent 
the approximate range over which each measurement applies. The full window functions in redshift space for each of these 
measurements are shown in panel (b); they have small leakage outside of the original redshift divisions. The window functions 
and the likelihoods are labeled in order of increasing redshift of the band powers in panel (a). The window functions satisfy 
three of our requirements: they make the band-powers uncorrelated, they are fairly well localized in redshift, and they are 
almost everywhere positive. In panel (a), we have used a uniform prior of 0.22 < Qm < 0.38 (a Gaussian prior Qm = 0.3 ±0.04 
gives very similar results), and we have assumed a flat universe throughout. 



dimensional parameter space (pi . . .J54, Qmi h). We al- 
low a generous range for the parameters pi (correspond- 
ing roughly to the vertical range in the left panels of the 
two Figures) and verify that changing the range leads to 
insignificant changes in the final constraints. We then 
marginalize the full likelihood over Qm and all values 
of h and project it onto the p\ . . .p4 space. We use a 
flat prior 0.22 < Om < 0.38, corresponding to the ±2cr 
allowed range from the joint analysis of various cosmo- 
logical probes Q. We have repeated our analysis with 
the Gaussian prior CIm = 0.30 ± 0.04 and found that 
the results are largely insensitive to the exact choice of 
either prior: the only notable change was that the first 
band power increases by about 0.15 with the Gaussian 
prior. The parameters pi are then rotated into the new 
parameters qi, which are now uncorrelated, following the 
methodology described in Sec. ITT1 

Figure 1 shows the final 68% and 95% CL constraints 
on the four band powers (i.e. the parameters qi) repre- 
senting w(z). We also show the weights that describe 
going from correlated parameters pi to the uncorrelated 
qi, as well as the full likelihoods of the four band pow- 
ers. The horizontal error bars in the left panel show 
the extent of the original bins; although the components' 
weights extend across the whole redshift range, the most 
weight (~ 60% or more) is in these respective bins and 
the band powers are therefore sufficiently localized in or- 
der to be easily interpreted. Note also that the weights 



are mostly positive and have small negative contribu- 
tions, as found in the context of matter power spectrum 
measurements |17| . 

As shown in Fig. 1, the equation of state is consistent 
with w = — 1 at the 95% CL in three out of four bins. 
We do find some (> 95% CL) evidence that w < — 1 at 
z < 0.2; however, to confirm this result with certainty 
will require more data, and in particular more stringent 
control of the systematic errors. Nevertheless, it is inter- 
esting that we find a similar tendency in the data as seen 
in completely independent analyses that use different, 
and less general, parameterizations |ll ll(|. The present 
approach, however, is less model dependent than these 
methods. In particular, any variations in the equation of 
state on redshift scales smoother than the binning scale 
can in principle be detected; more rapid oscillations can- 
not. This is why we consider this approach to be nearly 
model independent - it would be truly model indepen- 
dent if we used a large number of bins, as illustrated in 
Ref. [3. 

We now consider another parameterization of dark en- 
ergy - its energy density relative to the present value, 
f(z). We repeat the analysis and obtain constraints on 
f(z) shown in Fig. 2. They are roughly consistent with 
those for w, and are also consistent with the cosmolog- 
ical constant case at the 95% CL. Note that the likeli- 
hoods are fully contained in the allowed ranges, and we 
see no evidence for negative f(z). The weights of / are 
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FIG. 2: Same as Figure 1, but for f(z) = pde(z)/pde(0). Our band powers assume piecewise linear f(z), with the exception 
of the band power corresponding to the largest-redshift bin, which assumes constant f(z) across that bin. 



somewhat less well localized; however, the band powers 
are better determined than those of w, as expected from 
the fact that / is related to the luminosity distance data 
through a single, and not double, integral relation. 

IV. DISCUSSION AND SUMMARY 

We have used a variant of the principal component 
technique to produce uncorrelated, nearly model inde- 
pendent estimates of the equation of state of dark energy 
w(z) and its scaled energy density f{z). We used four 
redshift bins in each case, and found results that are in 
good agreement with previous analyses. We further ar- 
gued that the present approach nicely complements other 
methods that use conventional parameterizations of w(z) 
and f(z). Given that our band powers are uncorrelated, 
the interpretation of the cumulative evidence is particu- 
larly easy. 

If dark energy is due to the cosmological constant, then 
w = — 1 and f(z) = 1, and all of our band powers should 
be consistent with those values, independently of their 
window functions. Conversely, if we ever find strong sta- 
tistical evidence that even just one band power is differ- 
ent from —1 (for w) or 1 (for /), we will have ruled out 



the cosmological constant scenario. While we do find a 
hint of such evidence in the first band power of w{z), 
a definitive analysis will have to await more data and a 
careful assesssment of the systematics. 

While we have presented an analysis with ~ 150 su- 
pernovae and restricted ourselves to four bins in redshift, 
the generality and power of this method should make 
it perfectly suitable for the analysis of future supernova 
datasets, when the error bars are expected to improve 
by up to an order of magnitude and enable a much more 
quantitative analysis and comparison with models. Fur- 
thermore, the same techniques can be applied to a variety 
of other cosmological probes, as one can expect that their 
complementarity will considerably strengthen the SN la 
results. Finally, one can customize the proposed tech- 
nique specifically to maximize the return of any given 
test (say, whether w(z) = —1 or not). With an increase 
in the number of type la supernovae at high redshift, it 
is likely that these interesting possibilities will be consid- 
ered in the future. 
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